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Dependence of aj r c (inverse Sommer parameter in units of lattice spacing a) on am q (quark mass 
in lattice unit) has been observed in all lattice QCD simulations with sea quarks including the ones 
with improved actions. How much of this dependence is a scaling violation has remained an in- 
triguing question. Our approach has been to investigate the issue with an action with known lattice 
artifacts, i.e., the standard Wilson quark and gauge action with j3 = 5.6 and 2 degenerate flavors 
of sea quarks on 16 3 x 32 lattices. In order to study in detail the sea quark mass dependence, 
measurements are carried out at eight values of the PCAC quark mass values am q from about 
0.07 to below 0.015. Though scaling violations may indeed be present for relatively large am q , 
a consistent scenario at sufficiently small am q seems to emerge in the mass-independent scheme 
where for a fixed j3, 1/Vq and yfa have linear dependence on m q as physical effects similar to the 
quark mass dependence of the rho mass. We present evidence for this scenario and accordingly 
extract the lattice scale (a — 0.0805(7) fm, aT x = 2.45(2) GeV) by chiral extrapolation to the 
physical point. 
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1. Introduction 

The Sommer parameter, denoted by vq, is not a directly measurable quantity and as such may 
have uncertainties regarding its value. The method of determination of the lattice scale a using the 
Sommer parameter may also not be one's favorite way. However, the change of ro/a with am q , 
where m q is the sea quark mass, in simulations of all formulations (including improved versions) 
of lattice QCD has intrigued the lattice community for the last decade. Basically two questions 
appear: a) Is this a cut-off effect or a physical effect? b) How is the lattice scale a to be determined? 
Related issues are whether the scale a is to be taken as dependent on the quark mass m q . In that 
case, how does one chirally extrapolate hadronic quantities like masses given that quark masses are 
all at different scales? So far there is no theoretical understanding of the issues raised above. 



2. Our Simulation 



We have used unimproved Wilson gauge and fermion actions with well-known 0(a) cut-off 
effects with 2 degenerate light flavors of sea quark at a fixed gauge coupling given by j8 = 6/g 2 = 
5.6 on 16 3 32 lattices. One important feature of our investigation was large number, viz., 8 values 
of sea quark masses roughly in the range am q ~ 0.07 — 0.014. 

At each quark mass, 5000 equilibrated trajectories were generated using the standard HMC 
algorithm (DDHMC runs are underway on larger volumes). Using highly optimized gaussian 
smearing at both mesonic source and sink, pion and rho masses and their decay constants were 
computed [|]]. 

APE smearing was used on the gauge configurations with smearing level up to 40 with e = 2.5 
where c/4 = l/(e +4) was the coefficient of the staples. Expectation values of Wilson loops 
< W(R,T) > with rectangular extents R and T were then measured up to T = 16 and R = 8\/3. 
Reasonable plateau was obtained in effective potential versus T plots between T = 3 and T = 
5. The static potential aV(R) was extracted from single exponential fits between [T m [ ni T max ] = 
[3,4], [3,5], [4,5] using < W(R,T) >= C(R)exp[-aV(R)T]. Optimum smearing level was deter- 
mined at a given quark mass by observing the ground state overlap C(R) as a function of R (for 
details see [[T], §]). The optimum smearing level was found to be 30 for the lightest three quark 
masses and 25 for the rest. 

At each j8 and quark mass, the static potential aV (R) as obtained from the Wilson loops was 
analyzed with the phenomenologically successful well-known ansatz: 



aV(R) = aV Q + a 2 oR 



a 
R 



:OT 



where 5rot is the coefficient of the lattice correction term with 



4n ^ cos(aqfR) 
U q % 4sinW2) 



(2.1) 



(2.2) 



being the lattice fourier transform of the gluon propagator. 

The first 3 terms of aV (R) in eq. ^j] above is differentiated to obtain the Sommer parame- 
ter: a/r c = l/R c = ao l l 2 / \J (jV c — a) where .JVq = 1.65 and = 1 giving rise to the Sommer 
parameters ro and r\ respectively. 
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Figure 1: Fits of the static potential at K = 0.1575. The top panel shows the fitting range. 

3. Numerical Results 

A few general observations are noted regarding our fits of the static potential: (i) The difference 
([1/7?] — 1/7?) is never negligible on a finite lattice, (ii) a is expected to run with R at these 
intermediate length scales, (iii) We can only estimate an average a over the values of R where 
the static potential is fit. (iv) Perturbative running is generally applicable at scales > 2 GeV which 
translates into R < 1 in our case. 

We want to emphasize the importance of determining a, the coefficient of the 1 /R (Coulomb) 
term, origin of which is in the continuum perturbation theory. To determine it reliably, one naturally 
has to probe the small R region which has the problem of lack of rotational symmetry on the lattice. 
However, in our case, use of the correction term proportional to 5rot does the job as exemplified 
by fig. [I] where the fits describe the corrected data much beyond the fit range. However, to achieve 
such beautiful fits, one needs to tune all the fit parameters and the smearing level. Signature of 
good fits is not limited to fig. |l|. A good fit should also produce the values of 5rot close to that 
of a (unlike the random values as found in [J3j] ) and should show expected behavior of a, e.g., a 
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Figure 2: a versus R max plots at fixed values of R m i n and three T ranges for K = 0. 1575. 
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Figure 3: Left: a versus R m [ n at fixed 7? max for different values of the smearing level. Right: a versus am q 
at two values of R m ; n 



should increase as R m [ n increases and also as the smearing level increases (unlike wrong behavior 
as found in fig. 22 in [|[]). 

The above checks on the determination of a is crucial to our observation and inference. At 
Wilson hopping parameter K = 0.1575 {am q 0.03), Fig. || shows that a is relatively insensitive to 
change of /? max , however, the range [7 ^j n , T max ] = [3,4] and R m \ B = a/2 produces the most accurate 
determination of a. 

At the same K, the left panel of fig. || shows the behavior of a versus R m i n at various values 
of the smearing level and it exhibits the expected dependence on R m [ n and the smearing level. 
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Figure 4: ad 1 / 2 and a/r\ versus am q for two different values of R m i n 

It also shows that the most precise value of a is obtained at the smallest R m i n = a/2 and a is 
determined progressively imprecisely as R m i„ increases and beyond 7? m j n = 2 because of imprecise 
determination a does not grow with R m i n - 

The right panel of fig. || shows for a given R max = 6.708 and K = 0.1575 that a is weakly 
dependent on am q for two values of R m i n = y/2, 2. In particular, this dimensionless coefficient does 
not significantly depend on am q for small enough am q < 0.035. 

In fig. |] we show that for am q < 0.035, both aa 1 / 2 (left panel) and a/r\ (right panel) (and 
naturally also a/ro, although the plot is not shown here) can be fit linearly with am q : 

ao l/2 = C l +C 2 am q (3.1) 
a/r c = A c + B c am q (3.2) 

We also note that fig. f| shows data for two values of R m m used in previous figures. Qualitative 
conclusions about independence of a and linear dependence of ao 1 / 2 and a/r c on am q for am q < 
0.035 does not depend on the choice of R m { n as long as /? m ; n is small enough (obviously the data, 
especially a, is more accurate for smaller /? m i n )- In fact, the above qualitative conclusions do not 
depend on the choice of the values of the fitting parameters or the smearing level as long as they 
remain sensible. In addition, we observe from the right panel of fig. |]that although the individual 
values of a and ao 1 / 2 depend on the choice of R m i n , the value of a/r c at all am q is relatively 
insensitive to it. For our final analysis, we take [T m { n , T max ] = [3, 4], [R m i„, /? max ] = [V^, 3\/5] and 
APE smearing level is 30 for the lightest 3 quark masses and 25 for the rest of the quark masses. 
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Figure 5: Linear chiral extrapolation of a/r c and am p with (am n y 



Using the value of the lattice spacing (determined later) the smallest quark mass in use here is 
about 30 - 35 MeV and am q < 0.035 roughly translates into m q < 85 MeV at this j8 (= 5.6). Some of 
our smaller quark masses including the smallest one can be compared with values obtained with the 
same Wilson lattice QCD action parameters [§, |6|] but at larger lattice volumes and this comparison 
shows that our quark masses do not have any significant finite size (FS) effects. Our values of a/ro 
also are not expected to have any FS effect because the fitting range in R and the value r$ or r\ are 
well within the physical linear dimension of our lattice. 

Our quark masses are determined using PCAC and absence of FS effect on our quark masses 
establishes that PCAC is well satisfied on the lattice. The PCAC on the lattice differs from that in 
the continuum by 0(a). FS effect, if any, also enters through this 0(a) term. As a result, absence 
of FS effect is also an indirect indication of absence of scaling violations or cut-off effects. 



4. Interpretation of the Results 

Our numerical observation at fixed /3 of ara^-independence of the dimensionless coefficient 
a is interpreted as a signal for negligible cut-off effect in a for small enough am q . At fixed /3, 
the observed linear dependence of acr 1//2 and a/r c on am q for am q < 0.035 is then interpreted as a 
physical dependence of a 1 / 2 and \/r c onm q : 



a 1 / 2 
l/r c 



+ C 2 m q with C\ 
s^ c + B c m q with A c 



(4.1) 
(4.2) 
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In other words, for small enough am q , the cut-off dependence is negligible, a mass-independent 
scale-setting scheme at fixed /3 is implied, i.e., the lattice spacing a is fixed at a fixed /3, and the 
physical dependence of a 1 / 2 and l/r c on m q allows for a chiral extrapolation. 

Although absence of scaling violations cannot be solidly established unless one has data for 
different j8, the issue is that one still needs to set a scale at a fixed j3 and unless one has a criterion 
for getting rid of possible scaling violations, how is one ever going to achieve it and without that 
how would one do chiral extrapolation of hadronic observables? If any mass-independent scaling 
violations are present, that should then also invalidate all chiral extrapolations done to date on 
hadronic observables. One hopes, consistent with notions of universality, that at large enough j8 
(i.e., small enough lattice spacing a) and small enough quark mass m q , all valid formulations of 
lattice QCD should reach a regime where all cut-off effects are negligible. 

The left panel of fig. § then shows a linear chiral extrapolation of a/r c in (am„) 2 . The lat- 
tice spacing at the physical point is then extracted by solving a quadratic equation in the lattice 
spacing, along the lines of [§] (for details, see [0]). Our best determination of the lattice spacing 
is a = 0.0805(7) fm or a -1 = 2.45(2) GeV with r$ = 0.49 fm put in. This determination tallies 
very well with an independent hadronic determination of the lattice spacing through a linear chiral 
extrapolation of am p with {am n ) 2 (shown in the right panel of fig. |5|). The hadronic determination 
is less accurate and yields a = 0.0800(20) fm, i.e., cT x = 2.47(6) GeV. 

Details of this work can be found in [[7]]. 

Numerical calculations are carried out on a Cray XD1 (120 AMD Opteron@2.2GHz) sup- 
ported by the 10 th and 1 I th Five Year Plan Projects of the Theory Division, SINP under the DAE, 
Govt, of India. This work was in part based on the MILC collaboration's public lattice gauge 
theory code. See http://physics.utah.edu/ dtar/milc.html . 
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